各向异性两相材料尖劈奇性场的非协调元分析

提出了一个基于位移的、分析柱状各向异性两相材料尖劈端部邻域的奇性位移场和应力场问题的非协调元特征分析法. 该方法从柱状扇区的散度定理出发，将柱状扇区控制方程的弱式化为一个与虚功原理相同形式的方程，采用一种新的非协调元技术把所导出的``虚功原理'转化为标准一阶特征方程的求解问题. 非协调元法中，尖劈端部邻域的位移场假定没有采用奇异变换技术，有限元的单元形式是一维的. 将柱状各向异性两相材料尖劈视为``广义平面应变'问题，位移场与坐标z无关，只关注界面端的幂奇异性而不考虑对数奇异性. 运用该方法给出了柱状各向异性两相材料尖劈端部奇性应力指数、奇性位移角分布和应力角分布的算例. 所有的计算结果表明，该方法使用的单元少而且精度较高.     

Boundary contour analysis for surface stress recovery in 2-D elasticity and Stokes flow

Summary A variant of the boundary element method, called the boundary contour method, offers a further reduction in dimensionality. Consequently, boundary contour analysis of 2-D problems does not require any numerical integration at all. In a boundary contour analysis, boundary stresses can be accurately computed using the approach proposed in Ref. [1]. However, due to singularity, this approach can be used only to calculate boundary stresses at points that do not lie at an end of a boundary element. Herein, it is shown that a technique based on the displacement/velocity shape functions can overcome this drawback. Further, the approach is much simpler to apply, requires less computational effort, and provides competitive accuracy. Numerical solutions and convergence study for some well-known problems in linear elasticity and Stokes flow are presented to show the effectiveness of the proposed approach. This research was supported in part by the 2004 Ralph E. Powe Junior Faculty Enhancement Award from Oak Ridge Associated Universities and by the University of South Alabama Research Council.     

残余应力对压力容器破坏模式的影响分析

破前漏(简称LBB)是压力容器、核电站设备结构设计与评价中的一个重要准则.表面裂纹准静态扩展的几何形貌变化规律的预测是破前漏(LBB)评判十分重要的课题之一.本文对特定焊接残余应力场加载作用下,含三维表面裂纹的压力容器模型,用有限元软件(ABAQUS)进行了表面裂纹准静态扩展模拟计算,得到在此残余应力场作用下应力强度因子沿裂纹前缘的分布规律.结合外载引起的应力强度因子,就可以判别裂纹的扩展形貌,从而判断结构是否满足LBB要求.     

微分求积方法对于桩基大变形分析的应用

本文基于大变形的理论,采用弧坐标首先建立了具有初始位移的桩基的非线性数学模型,一组强非线性的微分-积分方程,其中,地基的抗力采用了Winkeler模型;其次,引入变数变换将微分-积分方程转化为一组非线性微分方程,并用微分求积方法离散了方程组,得到一组离散化的非线性代数方程;最后用Newton-Raphson迭代方法对离散化方程进行了求解,得到了桩基变形前后的构形、弯矩和剪力.计算中选取了两种不同类型的初始位移,并考察了它们对桩基大变形力学行为的影响.     

湍流最高激发态的局部结构及其绝对标度律

通过对湍流层次结构模型中提出的最高激发态的进一步研究发现，最高激发态 存在绝对标度律，且该绝对标度律是由信号中最强耗散涨落的局部结构产生的，并由此推测 出局部强间歇结构一般具有绝对标度的结论.     

Viscous waves and wave-structure interaction in a tank using adapting quadtree grids

Viscous waves and waves over a submerged cylinder in a stationary tank are simulated using a volume-of-fluid numerical scheme on adaptive hierarchical grids. A high resolution interface-capturing method is used to advect the free surface interface and the Navier–Stokes equations are discretised using finite volumes with collocated primitive variables and solved using a Pressure Implicit with Splitting of Operators (PISO) algorithm. The cylinder is modelled by using the technique of Cartesian cut cells. Results of flow of a single fluid past a cylinder at Reynolds number Re=100 are presented and found to agree well with experimental and other numerical data. Viscous free surface waves in a tank are simulated using uniform and quadtree grids for Reynolds numbers in the range from 2 to 2000, and the results compared against analytical solutions where available. The quadtree-based results are of the same accuracy as those on the equivalent uniform grids, and retain a sharp interface at the free surface while leading to considerable savings in both storage and CPU requirements. The nonlinearity in the wave is investigated for a selection of initial wave amplitudes. A submerged cylinder is positioned in the tank and its influence on the waves as well as the hydrodynamic loading on the cylinder is investigated.     

单向悬挂屋盖结构的风致气弹耦合效应数值模拟

综合运用计算流体力学和计算结构力学技术,建立了适用于索膜结构流固耦合风振分析的CFD数值模拟方法,并编制了相应的有限元计算程序。应用该程序对具有不同参数的大跨度单向屋盖结构进行了风振响应分析,探讨了来流风速、屋面质量和初始预张力等参数对结构流固耦合特性的影响。研究表明,由于屋盖前缘周期脱落的大尺度旋涡是诱导结构振动的主要原因,随着大尺度旋涡沿屋面向下游移动,其动能逐渐转化为结构变形能,使结构振动形态呈现明显的涡激振动特征。屋面质量、来流风速和初始预张力是影响结构流固耦合性能的主要参数。对于索膜结构,仅通过刚性模型来确定结构风压并进行风振响应分析,其结果是不可靠的。     

Micromechanics, macromechanics and constitutive modeling of the elasto-viscoplastic deformation of rubber-toughened glassy polymers

Under certain conditions, such as sufficiently low temperatures, high loading rates and/or highly triaxial stress states, glassy polymers display an unfavorable characteristic—brittleness. A technique used for reducing the brittleness (increasing the fracture toughness) of these materials is rubber toughening. While there is significant qualitative understanding of the mechanical behavior of rubber-toughened polymers, quantitative modeling tools for the large-strain deformation of rubber-toughened glassy polymers are largely lacking.In this paper, we develop a suite of numerical tools to investigate the mechanical behavior of rubber-toughened glassy polymers, with emphasis on rubber-toughened polycarbonate. The rubber particles are modeled as voids in view of their deformation-induced cavitation early during deformation. A three-dimensional micromechanical model of the heterogeneous microstructure is developed to study the effects of initial rubber particle (void) volume fraction on the underlying elasto-viscoplastic deformation mechanisms in the material, and how these mechanisms influence the macroscopic response of the material. A continuum-level constitutive model is developed for the large-strain elasto-viscoplastic deformation of porous glassy polymers, and it is calibrated against micromechanical modeling results for porous polycarbonate. The constitutive model can be used to study various boundary value problems involving rubber-toughened (porous) glassy polymers. As an example, the case of an axisymmetric notched bar is simulated for the case of polycarbonate with varying levels of initial porosity. The quality of the constitutive model calibration is assessed using a multi-scale modeling approach.     

三维电磁场辛有限元列式及电磁共振腔的计算

针对三维共振腔的电磁场分析,利用Maxwell方程的对偶方程体系形式,从其相应的对偶变量变分原理出发,导出了三维电磁场辛有限单元的详细列式。为了有限元列式的保辛,变分原理被积函数可导向对于对偶变量为对称的形式。变分原理的边界积分项对于相邻单元相互抵消。由于采用了对偶变量的插值函数,使得电磁场单元构造可以在层面上进行,从而避免了所谓的连续性问题。无物理意义的零本征解可采用奇异值分解加以排除。文末分别对矩形及圆柱形的共振腔做了数值计算并与解析解和棱边元计算结果进行对比,算例表明了列式及算法的有效性。     

点支式中空玻璃承载性能的有限元分析

中空玻璃承受均布横向荷载时,外、内片与密封气体之间产生复杂的交互作用,现行规范(程)基本按抗弯刚度比例在外、内片玻璃之间分配荷载。本文使用有限元法,分析点支式中空玻璃在承受均布横向荷载时的结构特性,利用迭代计算精确地求解外荷载在外、内片间的分配比例,计算与实验测值吻合得较好。本文指出,宽厚比、中空层厚度等因素对于荷载分配比例均有较大影响,在某些情况下,直接承受荷载的外片可能承受绝大部分的荷载,规范(程)的设计方法仍值得商榷。同时,本文建议在设计中考虑中空层温度变化对承载性能的影响。